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We calculate the s-wave scattering length and effective range and the p-wave scattering volume for 
^Li atoms interacting with atoms via the molecular potential. The length and volume 

are found by fitting the log-derivative of the zero energy wave function evaluated at short range 
to a long range expression that accounts for the leading van der Waals dispersion potential and 
then incorporating the remaining long range dispersion contributions to first order. The effective 
range is evaluated from a quadrature formula. The calculated parameters are checked from the zero 
energy limits of the scattering phase shifts. We comment on ill-conditioning in the calculated s-wave 
scattering length. 
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I. INTRODUCTION 

Du et alX presented the s-wave scattering length for 
^Li atoms interacting with ^33Cs atoms via the 
molecular potential. The elastic cross at very low ener¬ 
gies, needed to interpret the behaviour of an ultra-cold 
ensemble of atoms, is determined not only by the s-wave 
scattering length oq but also by the s-wave effective range 
Te and the p-wave scattering volume aii^ Scattering at 
higher angular momenta contributes little to the cross 
section. We calculated the length, range and volume. 

Du et alX used the interaction potential of Staanum et 
al^ We consider that generally such an interaction po¬ 
tential V (i?) consists of a short range part, consisting of 
a table of ab initio values at discrete values of the atomic 
separation R and an analytic expression that includes ex¬ 
change and dispersion, and a long range dispersion tail. 
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applicable from a separation beyond which exchange is 
negligible. In Eq. o Cn denotes a van der Waals coef¬ 
ficient; we regard only the first three dispersion terms as 
being important. 

With (fi denoting R times the radial parts of the zero- 
energy s-wave and p-wave functions with Z = 0,1 respec¬ 
tively (Z < 1 for this potential), we find 
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where /r is the reduced mass of the colliding atoms and 
h is Planck’s constant divided by 2tt. The asymptotic 
form of (j)i{R) depends on the length and the volume for 
Z = 0,1 respectively; with suitable normalisation 


The volume (Z = 1) is sometimes defined witout the 
factor three in Eq. We include the leading dispersion 
term —CqR~^ exactly and use the normalisation © to 
obtain a form of (j>i{R) that is valid at large Ri^ 
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where ag = 2717/ [r(l/4)]^, oi = 7^ [r(l/4)]^ /367r, T is 
the Gamma function, J, 21+1 are Bessel functions and 

^ 4 _ 

X = I2R^ with 7 = {/2\xCglW- . The quantities dg and 
cfi have dimensions of length and volume respectively and 
provide scaling for ag and oi; oi = 1.0642ao^ is very close 
to the cube of dg. The log-derivative of (/)j" (R) is 
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where the primes indicate differentiation of the Bessel 
functions with respect to the argument x. 

We calculated each scattering parameter ai as fol¬ 
lows. Using a recurrence relation^ we propagated the 
log-derivative ui{R) of the function (j)i{R) at short range 
out to a separation R* beyond which exchange is negli¬ 
gible; R* ~ 30 bohr in the present calculation. At long 
range we ignored the dispersion contribution —Cgi?”® — 
CigR~^'^ to the potential. By matching the calculated 
short range log-derivative ui{R*), which reflects all the 
potential, to the log-derivative uf {R) we obtained a value 
di{R*) of ai which is appropriate to the long range po¬ 
tential —CgR~^; the dependence of di{R*) on R* is a 
consequence only of the extent of the short range prop¬ 
agation with the complete potential. Having thus calcu¬ 
lated di{R*) we found o; from 
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where 
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accounts for the dispersion terms —CsR~^ and 
—CioR~^^. We calculated the Bessel functions and their 
derivatives by the method described by Press et ali^ We 
repeated the calculation with several values of R* to ob¬ 
tain convergence. This method allows calculations of oq 
and oi without the need to solve Eq. @ at large separa¬ 
tions. 

Equation © provides corrections to the length to 
first order in the dispersion contribution —C^R~^ and 
—with I — O.T’ By substituting Z = R^/3 when 
I = 1 we find that the volume is equivalent to the length 
for a new potential R~'^ x V{R) and a new zero-energy 
wave function R x (j)i{R)^ Therefore the volume cor¬ 
rection is determined by the length correction with the 
new potential. By substituting R~^ x V{R) for V{R), 
R X 4’i{R) for (j)o{R), Z for R^/3 and dZ for R^dR in 
Eq. ([7]) and then recasting the new quadrature in terms 
of R we obtain an equation for Si{R*) that is identical 
to Eq. © for / = 1. Hence Eq. © provides corrections 
for both length and volume. 

We found the range fromi^ 
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This quadrature based formula does not depend on the 
semiclassical approximation but follows the assumption 
that 1/R X (j)\{R) ^ good approximation to the zero- 

energy radial wave function throughout. This depends 
on two assumptions about the range: there is little in¬ 
fluence at small separations which was demonstrated 
semiclassically^ but can also be demonstrated quantum 
mechanically; and the dispersion terms —C^R~^ and 
—CiqR~^^ have little effect, which was demonstrated by 
perturbation calculationsj^ Thus the range for a pair of 
atoms depends almost solely on its scattering length oq 
and mean scattering length 

We, like Du et al.fi used the potential of Staanum et 
alji The values of ob and oi are 44.4028 bohr and 9.31722 
X10^ bohr^, respectively. We show the convergence of the 
length and volume in Table 1. With the approximatioit^ 
J^(x) ~ ^^j-Kx cos(x — i/7r/2 — 7r/4), we find: 
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and 


X*^ ^ X*^ , X* . , cos ^ — COSY 
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in which x* = 7^/2i?*^, 9 = (1 — |i|)i^7r, 0 = 2x* + ivii 
and X = with v = {21 + l)/4. Evaluation of ex¬ 
pressions (nni) and (HU is rapid but the values of the ap¬ 
proximate expression ([9]) converge less quickly than the 
accurate values as can be seen in Table 1; the departure 
from monotonic variation of the corrected lengths and 
volumes near 60 bohr is a consequence of the increasing 
inaccuracy of the trigonometric approximation with de¬ 
creasing X corresponding to increasing R. However the 
approximate expression provides a useful check on our 
results. 

Our final values of the length, volume and range are 
48.54 bohr, 1.617 xlO**’ bohr^ and 109.4 bohr, respec¬ 
tively. Our length differs by 4% from that, 50.5 bohr, 
calculated by Du et ah^ When the scattering length is 
large it is known that its calculated value can be sensitive 
to computational details and can even manifest spurious 
sign changes;^ but possible ill-conditioning in evaluations 
of various smaller scattering lengths is less well appreci¬ 
ated. Our reduced mass, 1.21481x10*'^ atomic units ob¬ 
tained from the atomic masses,— is essentially the same 
as that used by Du et alX but we found that our cal¬ 
culated ^Li-^^^Cs length varied by about 8% when the 
potential was changed by only 0.1%; a fractional reduc¬ 
tion of 0.05% (making the potential less negative beyond 
the classical turning point i?o) produced a length of 50.5 
bohr. Thus it is plausible that our result and that of Du 
et ali^ differ because of possible different interpolations 
used to calculate the potential. With the nuclear masses 
we found the scattering length to be 49.44 bohr. 

The cause of the sensitivity of the calculated length 
can be seen in the semiclassical expression for it^ 


Og® = og [1 — tan(<i) — 7r/8)] (12) 

where fi<l> is the action integral: 
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The angle $ — n/S has value 172.702. The potential 
supports 55 bound states and the argument of the tan¬ 
gent function in Eq. m is the same as $ — 7r/8 — BBtt = 
—0.085. A 0.1% change in the potential alters the action 
integral by 0.05% and changes $ by 0.86 which alters 
$ — 7r/8 — 557r by 10%. The semiclassical length is 48.20; 
the value of $ — tt/S — 557r that allows Eq. ((T7|) to repro¬ 
duce the quantal length is 0.093 for which $ — 7r/8 — SStt 
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TABLE I. Scattering lengths (bohr) and volnmes (10“^ bohr®) 
for cold ^Li - ^®®Cs collisions, for varions valnes of separation 
R* (bohr). 


R* 

do 

ucu 


a'l 

“U 

aj 


30 

53.68 

48.36 

48.28 

4.642 

4.921 

1.030 

1.122 

35 

52.85 

48.37 

48.31 

4.098 

4.379 

1.199 

1.269 

40 

50.68 

48.48 

48.37 

2.883 

3.069 

1.493 

1.601 

45 

49.25 

48.54 

48.48 

2.056 

2.281 

1.602 

1.664 

50 

48.78 

48.54 

48.53 

1.757 

2.037 

1.615 

1.632 

55 

48.68 

48.54 

48.55 

1.675 

1.987 

1.617 

1.611 

60 

48.67 

48.54 

48.56 

1.657 

1.981 

1.617 

1.601 

65 

48.67 

48.54 

48.56 

1.655 

1.981 

1.617 

1.599 

100 

48.59 

48.54 

48.54 

1.641 

1.941 

1.617 

1.614 

500 

48.54 

47.54 

48.54 

1.617 

1.915 

1.617 

1.617 


^ With corrections I©. 0. 

^ With approximate corrections ®, 0, Gnii, ETj. 
^ Predicted by Gao’s formula flit . 


changes by 9% which is consistent with our numerical 
prediction. 

We made a simple check on our calculations by com¬ 
paring the values of di with those predicted from do by 
Gao’s formulapi^ 


0-0 ~ 0,0 
ai = ai X ^ 

oo - ao/2 


(14) 


which applies to potentials as i? ® as i? —> oo at 
separations where the Bessel functions of Eq. are 
well represented by the trigonometric approximation and 
the rotational term 2/i?*^ is small. Table 1 shows good 
agreement between the volumes predicted by Eq. o 
for i? Ri 30 bohr, where these conditions are met, but the 
agreement deteriorates for larger values of R as expected. 

As a check on all our calculated parameters we solved 
the s-wave and p-wave scattering equations over a large 
but finite range of R at low energies and fitted the scat¬ 
tering parameters to the effective range expansions of the 
phase shifts^ We extrapolated the results to zero energy 
but in doing so we included only even powers of the wave 
number which is legitimate for a truncated potential but 
the long range nature of the actual potential must be ac¬ 
counted for.— We carried out this accounting by adding 
corrections4dd We found the same values for the volume 
and length. The range thus found accounts exactly for 
the dispersion terms —Cgi?”® and —CioR~^^; our value 
changed only slightly to 109.6 bohr, verifying that these 
terms have little influence.— 

In conclusion, we have calculated the scattering param¬ 
eters of length, volume and range for ^Li atoms interact¬ 
ing with ^®®Cs atoms via the molecular potential 

as 48.54 bohr, 1.617 xlO"^ bohr® and 109.4 bohr, respec¬ 
tively. 
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